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When the universe was about 10 /^seconds old, a first order cosmological quark - hadron phase 

transition occurred at a critical temperature of around 200 MeV. In this work, we study the quark- 

hadron phase transition in the context of brane-world cosmologies, in which our Universe is a 

three-brane embedded in a five-dimensional bulk, and within an effective model of QCD. We an- 

QO , alyze the evolution of the physical quantities, relevant for the physical description of the early 

(^ • universe, namely, the energy density, temperature and scale factor, before, during, and after the 

("^ ' phase transition. To study the cosmological dynamics and evolution we use both analytical and 

OQ . numerical methods. In particular, due to the high energy density in the early Universe, we consider 

in detail the specific brane world model case of neglecting the terms linearly proportional to the 

energy density with respect to the quadratic terms. A small brane tension and a high value of 

the dark radiation term tend to decrease the effective temperature of the quark-gluon plasma and 

^\ ■ of the hadronic fluid, respectively, and to significantly accelerate the transition to a pure hadronic 

phase. By assuming that the phase transition may be described by an effective nucleation theory, 

we also consider the case where the Universe evolved through a mixed phase with a small initial 

supercooling and monotonically growing hadronic bubbles. 
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Kj^ . I. INTRODUCTION 

^. 

^^^ • The possibihty that our 4D universe may be viewed as a brane hypersurface embedded in a higher dimensional 

^—^ ' bulk space has attracted considerable interest lately. A scenario with an infinite fifth dimension in the presence of 
. . a brane can generate a theory of gravity which mimics purely four-dimensional gravity, both with respect to the 
C ■ ' classical gravitational potential and with respect to gravitational radiation [1]. The gravitational self-couplings are 
^^ , not significantly modified in this model. This result has been obtained from the study of a single 3-brane embedded 
~!~ in five dimensions, with the 51? metric given by ds^ = e~^^^'r]^iydx^dx^ -f dy^ , which can produce a large hierarchy 

T\ , between the scale of particle physics and gravity due to the appearance of the warp factor [1]. Even if the fifth 
^ • dimension is uncompactified, standard AD gravity is reproduced on the brane. In contrast to the compactified case, 
this follows because the near-brane geometry traps the massless graviton. Hence this model allows the presence of 
large or even infinite non-compact extra dimensions. Our brane is identified to a domain wall in a 5-dimensional 
cd ' anti-de Sitter space-time. 

The effective gravitational field equations on the brane-world, in which all the matter forces except gravity are 
confined on the 3-brane in a 5-dimensional space-time with Z2-symnietry have been obtained, by using an elegant 
geometric approach, in [2, 3]. The correct signature for gravity is provided by the brane with positive tension. If 
the bulk space-time is exactly anti-de Sitter, then generically the matter on the brane is required to be spatially 
homogeneous. The electric part of the 5-dimensional Weyl tensor Ejj gives the leading order corrections to the 
conventional Einstein equations on the brane. This implies a modification of the basic equations describing the 
cosmological and astrophysical dynamics, which has been extensively considered [4]. For reviews of the dynamics and 
geometry of the brane Universes, as well as for the discussions of the cosmological implications see [5]. 

According to standard cosmology, as it expanded and cooled, the early Universe is expected to have undergone 
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a series of symmetry-breaking phase transitions, at which topological defects may have formed. Phase transitions 
are labelled first or second order, according to whether the position of the vacuum state in field space changes 
discontinuously or continuously, as the critical temperature is crossed. A first order phase transition proceeds by 
bubble nucleation and expansion. When at least (4 — n) of these bubble collide, for n — 0, 1,2, an n-dimensional 
topological defect may form in the region between them [6]. Recent lattice QCD calculations for two quark flavors 
suggest that QCD makes a transition at a temperature of Tc ^ 150 McV [7]. This phase transition, which may have 
occurred in the early Universe, could lead to the formation of relic quark-gluon plasma objects, which still survive 
today. 

A first order quark- hadron phase transition in the expanding Universe can be described generically as follows [6] . As 
the color deconfined quark-gluon plasma cools below the critical temperature T^ « 150 MeV, it becomes energetically 
favorable to form color confined hadrons (primarily pions and a tiny amount of neutrons and protons, due to the 
conserved net baryon number). However, the new phase does not show up immediately. As is characteristic for a 
first order phase transition, some supercooling is needed to overcome the energy expense of forming the surface of 
the bubble and the new hadron phase. When a hadron bubble is nucleated, latent heat is released, and a spherical 
shock wave expands into the surrounding supercooled quark-gluon plasma. This reheats the plasma to the critical 
temperature, preventing further nucleation in a region passed by one or more shock fronts. Generally, bubble growth 
is described by deflagrations, with a shock front preceding the actual transition front. The nucleation stops when 
the whole Universe has reheated to Tc- This part of the phase transition passes very fast, in about 0.05 /xs, during 
which the cosmic expansion is totally negligible. After that, the hadron bubbles grow at the expense of the quark 
phase and eventually percolate or coalesce. The transition ends when all quark-gluon plasma has been converted to 
hadrons, neglecting possible quark nugget production. The physics of the quark-hadron phase transition, as well as 
the cosmological implications of this process have been extensively discussed in the framework of general relativistic 
cosmology in [8] - [21] 

In the context of brane- world scenarios, the Friedmann equation contains deviations to the 4Z? case, which results 
in an increased expansion in early times. This in general has important effects on the cosmological evolution, and in 
particular on cosmological phase transitions. First order phase transitions have also been considered, in the framework 
of the brane- world model, in [22]. Due to the effects coming from the extra-dimensions, phase transitions require a 
higher nucleation rate to complete, and baryogenesis and particle abundances could be suppressed. The evolution of 
topological defects is also affected, but the increased expansion cannot solve the monopole and domain wall problems. 

In this work, we consider the quark-hadron phase transition in the brane-world scenario. By assuming that the 
phase transition is of the first order, we study in detail the evolution of the relevant cosmological parameters (energy 
density, temperature, scale factor, etc) of the quark-gluon and hadron phases, and the phase transition itself. It is 
important to emphasize that in the early universe the energy density is high, so that one may neglect the terms 
linearly proportional to the energy density with respect to the quadratic terms. This is carried out in detail in this 
work. An important parameter to describe the phase transition is the hadron fraction, whose time evolution describes 
the conversion process. We also consider the effect of the dark radiation term on the phase transition. The brane 
world effects (the quadratic corrections to the matter energy-momentum tensor, described by the numerical value 
of the brane tension, and the dark radiation term) lead to an overall decrease of the temperature of the very early 
universe, and accelerate the transition to the pure hadronic phase. 

This paper is organized in the following manner. In Section II, we briefly outline, for self-completeness and self- 
consistency, the field equations in brane-world models and the basics of the quark-hadron phase transition. In the 
latter, we lay down the equations of state and the relevant physical quantities that are analyzed in the remaining 
Sections. In Section III we analyze in detail the quark-hadron phase transition. In Section IV, we consider bubble 
nucleation in the brane-world scenario, by assuming that the phase transition may be described by an effective 
nucleation theory. We discuss and summarize our results in Section V. In the present paper we use a system of units 
so that the speed of light is c = 1. 

II. GEOMETRY, BRANE-WORLD FIELD EQUATIONS, EQUATIONS OF STATE AND 

CONSEQUENCES 

A. The field equations in the brane-world models 

In the brane-world models the effective four-dimensional gravitational field equations on the brane take the form 

[2,3]: 



where A = fc|(A 
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T^i, is the matter energy-momcntuni tensor on the brane and 



T = T^ ^ is the trace of the energy-momentum tensor. The Einstein equation in the bulk also implies the conservation 
of the energy momentum tensor of the matter on the brane [3]. The first correction term relative to Einstein's general 
relativity is the inclusion of a quadratic term 5*^,^ in the energy-momentum tensor, arising from the extrinsic curvature 
terms in the projected Einstein tensor, and is given by 
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The second correction term, E^^, is the projection of the 5-dimensional Weyl tensor, C'abcDj onto the brane, and is 
defined as E"^^ ~ ^u^v ^ABCD'^^n^i and encompasses nonlocal bulk effects. The only general known property of 
this nonlocal term is that it is traceless, i.e., -E^^ = 0. The symmetry properties of iJ^^ imply that in general we can 
decompose it irreducibly with respect to a chosen 4-vclocity field u^ as 
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where h^i, = .g^/y + u^u^ projects orthogonal to u^, U is a scalar, Q^ a spatial vector and Vf^i, a spatial, symmetric and 
trace-free tensor, respectively. For homogeneous and isotropic Friedmann-Robertson- Walker (FRW) type cosmological 
models Q^ = Vf^u = [5], and hence the only non-zero contribution from the 5-dimensional Weyl tensor from the 
bulk is given by the scalar term U. 

There are several constraints which have been obtained for the brane tension A. Thus from big bang nucleosynthesis 
constraints it follows A > 1 MeV* [23] . A much stronger bound for A arises from null results of submillimcter tests of 
Newton's law, giving A > 10^ GeV^ [24]. An astrophysical lower limit on A, which is independent of the Newton law 
and the cosmological limits has been derived in [23], leading to A > 5 x 10^ MeV^. 

In this work, we assume that the space-time geometry is the flat FRW metric, given by 



ds^ = -dt^ + a'^{t) (dx^ + dy^ + dz^) 



(4) 



For the matter energy-momentum tensor on the brane we restrict our analysis to the case of the perfect fluid energy- 
momentum tensor. 



T^" = {p + v)u''u'' + pg""" , 

where p and p are the energy density and isotropic pressure of the cosmological fluid on the brane, respectively. 
Thus, the gravitational field equations, corresponding to the line element (4) become 
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Eq. (9) can be immediately integrated to give the following general expression for the "dark energy" U: 
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with Wo > a constant of integration. 



B. Quark- hadron phase transition 



In this Section, we outline the relevant physical quantities of the quark-hadron phase transition, which will be 
applied in the following sections in the context of the brane-world scenario. Note that the scale of the cosmological 
QCD transition is given by the Hubble radius R[j at the transition, which is Rh ~ "m-pi/T^ ~ 10km, where T^. is 
the critical temperature. The mass inside the Hubble volume is ^ IMq. The expansion time scale is 10~^ s, which 



should be compared with the timescale of QCD, 1 fm/c« 10~^^ s. Even the rate of the weak interactions exceeds the 
Hubble rate by a factor of 10^. Therefore, in this phase photons, leptons, quarks and gluons (or pions) are lightly 
coupled and may be described as a single, adiabatically expanding fluid [12]. 

In order to study the quark-hadron phase transition it is necessary to specify the equation of state of the matter, 
in both quark and hadron state. Giving an equation of state is equivalent to give the pressure as a function of the 
temperature T and chemical potential /i. At high temperatures the quark chemical potentials are equal, because weak 
interactions keep them in chemical equilibrium, and the chemical potentials for leptons are assumed to vanish. Thus 
the chemical potential for a baryon is defined by /is = 3/ig. The baryon number density of an ideal Fermi gas of three 
quark flavors is given by ns ~ T^/xb/3, leading to hb/T ~ 10~^ at T > Tc- At low temperatures hb/T ^ 10^^. 
Therefore the assumption of a vanishing chemical potential at the phase transition temperature in both quark and 
hadron phase represents an excellent approximation for the study of the equation of state of the cosmological matter 
in the early Universe. In addition to the strongly interacting matter, we assume that in each phase there are present 
leptons and relativistic photons, satisfying equations of state similar to that of hadronic matter [6] . 

The equation of state of the matter in the quark phase can generally be given in the form 

p,^3agT^ + V{T), p,=aqT^-V{T), (11) 

where a^ = (tt^/QO) .g,, with .g, = 16 + (21/2)iVi^ + 14.25 = 51.25 and Np = 2. V (T) is the self-interaction potential. 
For V we adopt the expression [19] 

V{T) = B + jtT^ - utT^, (12) 

where B is the bag constant, ax ~ 77r^/20, and "fr — Ws/4, with TOs the mass of the strange quark in the range 
rus € (60 — 200) MeV. The form of the potential V corresponds to a physical model in which the quark fields are 
interacting with a chiral field formed with the tt meson field and a scalar field. If the temperature effects can be ignored, 
the equation of state in the quark phase takes the form of the MIT bag model equation of state, Pq = {pq — 4i?)/3. 
Results obtained in low energy hadron spectroscopy, heavy ion collisions and phcnomcnological fits of light hadron 
properties give B^/^ between 100 and 200 MeV [25]. 

In the hadron phase we take the cosmological fluid to consist of an ideal gas of massless pions and of nuclcons 
described by the Maxwell-Boltzmann statistics, with energy density ph and pressure pt, respectively. The equation 
of state can be approximated by 

Ph (T) = ^Ph (T) = a^T\ (13) 

where a^r = (tt^/QO) gh and gh = 17.25. 

The critical temperature Tc is defined by the condition pq (Tc) = Ph [Tc) [6], and is given, in the present model, by 



rj. ^ I IT + Vtt + 4 (Qg + ar - a^) B 
'V 2iaq + aT-a^) ' ^ ' 

For rris = 200 MeV and B^/'^ = 200 MeV the transition temperature is of the order T^ w 125 MeV. Since 
the phase transition is of first order, all the physical quantities, like the energy density, pressure and entropy exhibit 
discontinuities across the critical curve. The ratio of the quark and hadron energy densities at the critical temperature, 
Pq {Tc) IPh {Tc), is of the order of 3.62 for m^ = 200 MeV and B^^^ = 200 MeV. If the temperature effects in the 
self-interaction potential V are neglected, ax = 7t ~ 0, then we obtain the well-known relation between the critical 
temperature and the bag constant, B ~ {gq — gh) 7r^T^/90 [6]. 

III. DYNAMICS OF THE BRANE UNIVERSE DURING THE QUARK-HADRON PHASE 

TRANSITION 

The quantities to be traced through the quark-hadron phase transition in the brane-world cosmological scenario are 
the energy density p, the temperature T and the scale factor a. These quantities are determined by the gravitational 
field equations (6) and (8) and by the equations of state (11), (12) and (13). We shall consider now the evolution of 
the brane-world before, during and after the phase transition. 

Before the phase transition, T > Tc the brane-world is in the quark phase. With the use of the equations of state 
of the quark matter, and the Bianchi identity on the brane, Eq. (8) can be written in the form 
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and can be integrated to give the following scale factor-tcnipcrature relation: 
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where gq is a constant of integration. 

With the use of Eq. (15), from the gravitational field equations we obtain an expression, describing the evolution 
of the temperature of the branc Universe in the quark phase, given by 



dT 
'dt 



T3 JAoT» + AiT^ + ^2^4 + A3T2 + AiT^A^ exp (-2^6^) + A-j 
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V3 A^T^ + Aq 

where we have denoted 

Ao = kl (3a, - arf /2A, Ai ^ kj (Sa, - ar) 7t/A, A2 = kj [-f^/2X + (Sa, - ar) (1 + B/X)] , (18) 
^3-fc|(l + S/A)7T, Ai^6Uo/\kl4,A5:^i3ag-aT)/3ag, Ae=-fT/6ag, A7 = A + fc|B (1 + B/2A) . (19) 

The variation of the temperature as a function of the parameter t = k/^t in the quark matter filled branc-world is 
represented, in the case of a vanishing cosmological constant A = 0, for different values of the brane tension A and for 
a fixed numerical value of the constant A4 = 6ZYo/Afc|aQ, in Fig. 1. 
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FIG. 1: Variation of the temperature of the quark fluid on the brane as a function of r = k4t for different values of the brane 
tension A; A = 5 x 10** MeV* (solid curve), A = 5 x 10^ MeV* (dotted curve), A = 5 x 10^° MeV* (short dashed curve) and 
A = 5 X 10^'' MeV* (long dashed curve). We have assumed that in this phase the cosmological constant is vanishingly small. 
For the bag constant we have chosen the value B^'** — 200 MeV, while A4 — 0.01. 



The variation of the temperature of the brane in the quark phases for different values of the constant A4, corre- 
sponding to different numerical values of the integration constant Uq, describing the effect of the dark radiation term 
on the cosmological evolution, and considering a fixed value for the brane tension A, is represented in Fig. 2. 

In order to have an analytical insight into the evolution of the cosmological quark matter in the brane-world and 
on the effects of the extra-dimensions on the phase transition, we consider the simple case in which the temperature 
corrections can be neglected in the self-interaction potential V. In this case V = B = constant, and the equation 
of state of the quark matter is given by the bag model equation of state, pq = (pg — 4B) /3. Thus, Eq. (8) can 
immediately be integrated to give the following scale factor-temperature relation: 



a{T) = — , 



(20) 



where qq is a constant of integration. Since in the initial stages of the evolution of the Universe the density of the 
matter is very high, we suppose that in Eq. (6) we can neglect the term proportional to the energy density with 
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FIG. 2: Variation of the temperature of the quark fluid on the brane as a function of r = kit for different values of the dark 
radiation term U: A4 = 10^ (solid curve), A4 = 10* (dotted curve), A4 = 10* (short dashed curve) and A4 = 10® (long dashed 
curve). We have assumed that in this phase the cosmological constant is vanishingly small. For the bag constant and for the 
brane tension we have chosen the values B^''*' = 200 MeV and A = 5 x 10® MeV*, respectively. 



respect to the term containing p^. The contributions of the cosmological constant and of the dark energy are also 
negligible small. Therefore the evolution of the quark phase of the brane- world is described by the equation 
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Hence the time dependence of the temperature can be obtained from the equation 
vi^ith the general solution given by 



(21) 



(22) 



T{t) 



51/4^1/4 exp(-^f 



1 — 3a„Cexp 



4^:45 
' Vex 



1/4 ■ 



(23) 



The integration constant C is related to the temperature Tq of the quark matter at the time tg by means of the 
relation 
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leading to 
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During the phase transition, the temperature and the pressure are constants, T — T^ and p — pc, respectively. 
The entropy 5* = sa^ and the enthalpy W = {p + p) a^ are conserved quantities. During the phase transition p (t) 
decreases from pq (T^) = pg to pi^ (T^) = pn. For phase transition temperature of T^ = 125 MeV we have pg « 5 x 10^ 
MeV' and p^ w 1.38 x 10^ MeV'', respectively. For the same value of the temperature the value of the pressure of 



the cosmological fluid during the phase transition is pc w 4.6 x 10* MeV*. It is convenient, following [6], to replace 
p{t) by h{t), the volume fraction of matter in the hadron phase, by defining 



p (t) = pHh{t) + PQ [1 - h{t)] = PQ [1 + nh{t)] , 



(26) 



where we denoted n — [pn ~ Pq) / Pq- At the beginning of the phase transition h{tc) = 0, where tc is the time 
corresponding to the beginning of the phase transition, and p {tc) = PQ, while at the end of the transition h [th) = 1, 
where f^ is the time at which the phase transition ends, corresponding to p{th) = Ph- For t > tf^ the Universe enters 
in the hadronic phase. 
From Eq. (8) we obtain 
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where we denoted r = [pn — pq) / {pq + pc)- The above equation immediately leads to 
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where we have used the initial condition h (tc) = 0. The evolution of the fraction of the matter in the hadron phase 
is described by the equation 
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It is also convenient to plot the variation of the hadron fraction as a function of the different parameters. The 
variation of the hadron fraction /i as a function of the parameter r — k^t is represented, for different values of the 
brane tension and for a fixed value ofUo in Fig. 3. 
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FIG. 3: Variation of the hadron fraction h as a, function of the parameter r = k4t during the quark-hadron phase transition 
on the brane, for different values of the brane tension: A = 5 x 10* MeV* (solid curve), A = 5 x 10^ MeV* (dotted curve), 
A = 5 X 10^" MeV* (short dashed curve) and A = 5 x 10^'^ MeV* (long dashed curve). We have assumed that the effect of the 
cosmological constant can be negfected. The temperature during the phase transition is Tc — 125 MeV. For the bag constant 
we have chosen the value B^^* = 200 MeV, while ^4 = 0.01. 



The time variation of the hadron fraction during the phase transition for a fixed brane tension A and for different 
values of the integration constant Uq is represented in Fig. 4. 

If the quadratic term in the energy density, describing the effects of the extra-dimensions, dominates the evolution 
of the Universe, and by neglecting the dark energy contribution, the evolution of the hadron fraction during the phase 
transition is described by the equation 
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FIG. 4: Variation of the hadron fraction h as & function of the parameter r = k4t during the quark-hadron phase transition 
on the brane, for different values of the dark radiation term uo = SUo/K^Xa* (tc): Mo ~ 10^'^ MeV^ (sohd curve), uo = 10^^ 
MeV^ (dotted curve), mo ~ 10^^ MeV^ (short dashed curve) and mo — 10^^ MeV^ (long dashed curve). We have assumed that 
the effect of the cosmological constant can be neglected. The temperature during the phase transition is Tc — 125 MeV. For 
the bag constant we have chosen the value B^''' = 200 MeV, while A = 5 x 10*. 



with the general solution given by 
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where we have also used the initial condition ft. (tc) — 0. 

The respective evolution of the scale factor during the phase transition for the extra-dimensions dominated evolution 
of the cosmological fluid is described by the equation 
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The phase transition ends when h{t) has increased to 1. If the evolution is dominated by the terms coming from 
the extra-dimensions in the bulk, the time th at which the phase transition ends is 
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At the end of the phase transition the scale factor of the Universe has the value a (th) = a (tc) (r + 1) ' . 

Finally, after the phase transition, the energy density of the pure hadronic matter is ph = iph = iciTrT'^- The 
Bianchi identity Eq. (8) gives a{T) = a (th) Tc/T. The temperature dependence of the brane Universe in the hadronic 
phase is governed by the equation 
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As before, it is convenient to plot the variation of the temperature in terms of the various parameters. The variation 
of the temperature of the hadron fluid filled brane Universe as a function of the brane tension A is represented in Fig. 
5. 

The time variation of the temperature for the brane in the hadron phase for different values of the dark energy 
coefficient Uq is represented in Fig. 6. 
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FIG. 5: Variation of the temperature of the hadron fiuid on the brane, as a function of r = k4t, for different values of the 
brane tension A; A = 5 x 10* MeV* (solid curve), A = 5 x lO'' MeV* (dotted curve), A = 5 x 10^° MeV* (short dashed curve) 
and A = 5 X 10^^ MeV* (long dashed curve). We have assumed that in this phase the cosmological constant is vanishingly 
smaU. For the bag constant we have chosen the value B^''' = 200 MeV, while the dark radiation term has been fixed so that 
6Uo/\kia^ (th) T^ = 0.01. 
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FIG. 6: Variation of the temperature of the hadronic fluid on the brane, as a function of r = kit, for different values of the 



dark radiation term uo = GUo/Xk^a (th)T^ ■ Mo = 10 (solid curve), uo = 10 (dotted curve). 



Uo 



10 (short dashed curve) 



and Mo ~ 10^^ (long dashed curve). We have assumed that in this phase the cosmological constant is vanishingly small. For 
the bag constant and for the brane tension we have chosen the values B^'* = 200 MeV and A = 5 x 10* MeV*, respectively. 



If the extra-dimensional terms dominates the four-dimensional gravitational effects, then the time variation of the 
temperature of the cosmological hadron fluid is given by 



T{t) 



12k^a-^ 



-1-1/4 



6A 



{t - th) + r; 



(35) 



where we have used the initial condition T {th) = T^. 

In the extra-dimensions dominated hadron phase the scale factor of the Universe evolves according to the law 



a{t) Ri a (th) Tc 



12kiia-^ 



6A 



{t-th) + T- 



1/4 



(36) 
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IV. BUBBLE NUCLEATION IN THE BRANE- WORLD COSMOLOGICAL SCENARIOS 



In order to describe the process of formation and evolution of microscopic quark nuclei in the cosmological fluid on 
the brane one must use nucleation theory. The goal of nucleation theory is to compute the probability that a bubble 
or droplet of the A phase appears in a system initially in the B phase near the critical temperature [26] . Homogeneous 
nucleation theory applies when the system is pure. Nucleation theory is applicable for first-order phase transitions 
when the matter is not dramatically supercooled or superheated. If substantial supercooling or superheating is present, 
or if the phase transition is second order, then the relevant dynamics is spinodal decomposition [27]. 

In this context, note that in the analysis of the phase transition considered in the previous section we do not take 
into account the possibility that the quark plasma could undergo a supercooling phase, in which it stays in the quark 
phase below the critical temperature. If the supercooling is not too strong, the phase transition can be described by 
an effective nucleation theory, as mentioned above. As the temperature decreases, there is a probability that a droplet 
of hadrons is nucleated from the quark plasma. The nucleation probability density is given by [6] 



P{t) = PoTc cxp 



Wq 



(l-T^y 



(37) 



where T — T/Tc- Once a droplet of hadrons is formed, it starts to expand explosively [28], with a velocity Vfr smaller 
than the sound speed. Contemporarily, a much quicker shock wave is generated. More and more bubbles are created 
while the temperature decreases, until the shock waves collide and re-heat the plasma to the critical temperature. 

To calculate the fraction of the volume which at a time t is turned to the hadronic phase in the small supercooling 
scenario, we should sum over the volumes at the time t of the bubbles which are nucleated at a previous time tp, 
times the probability to create a bubble at that time tp, which is given by 



fit) 



dtpPWYb^frit-tpn 



(38) 



where U is the time at which the critical temperature is reached, and the nucleation process started. 
(22)-(23), we can express this integral in terms of normalized temperature. After some algebra we get 



Using Eqs. 



fiT) 



dT, 






6A 



2 exp 



i^-T^r 



Tp (3a,r4 + b) 



log 



'J^q-'-p 



m^ 



"^^q^p 



(39) 



where b = B/T^. 

The integrand is extremely peaked around its maximum, so the integral is dominated by the value of the function 
at the maximum, which can be found by solving the equation: 



log 



iaqTp 



b{^ 



'ia^T^ 



SwqT^ 
(l-Tp^y 



7^2 

p 



3agT^ + I 




3agT4 -I- I 



(40) 



To evaluate this maximum, we note that the integrated function falls very rapidly to zero as Tj, — > 1 because of the 
exponential, so the maximum must be very close to the lower end of the interval. Thus we can set Tp^max ~ T -I- a; 
with x/T <C 1 and approximate the first derivative (40) to first order in x. Moreover, since we will assume only a 
small supercooling, the final temperature T is very close to 1, so that the first term in the square bracket of the left 
hand side in Eq. (40) is much larger than the other two, which arc of order 1 and can be discarded. Eventually we 
obtain 



x{f) 



(1 - f*Y {6agf* + b) 
8wo6T3 



(41) 



so that the integral (39) can be approximated as 



/(t) c, x{f) 



Povj^ 



6A 



exp 



\l-{T+x{T)y 



B^ V«l/ if + xif))[3aq(f + x{f))^ + b] 



log 



3ag{T + x{T)y 



T+x{T) 
T 



3ag{T + x{T)y + b 



(42) 
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FIG. 7: Variation of the hadronic fraction / as a function of the normahzed temperature in the small supercooling scenario, 



for different values of the brane tension A: A = 5 x 10* MeV* (solid curve), A 
MeV* (short dashed curve) and A = 5 x 10^'' MeV* (long dashed curve). The t 
the values: B'^^* = 200 MeV and Tc = 125 MeV, respectively. 
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FIG. 8: Variation of the hadronic fraction / as a function of the parameter r = kit in the small supercooling scenario, for 
different values of the brane tension: A = 5 x 10* MeV* (solid curve), A = 5 x 10^ MeV" (dotted curve), A = 5 x 10^° MeV* 
(short dashed curve) and A = 5 x 10^^ MeV* (long dashed curve). The bag constant and the critical temperature have the 
values: B^/* = 200 MeV and Tc = 125 MeV, respectively. 



The plot of /(T) is presented in Fig. 7 for values of the parameters used previously, and setting the front velocity 
Vfr — 10^^ and the other constants Wo,Po = 1- The fraction of hadronic matter stays very close to zero until the 
supercooling temperature is between T — 0.98 and T — 0.97, then it jumps to 1 very rapidly. The same behavior can 
be verified in the time dependent plot of Fig. 8. The comparison between Fig. 8 and Fig. 3 shows clearly the difference 
between the first order phase transition and the supercooling. In the former, the hadronic fraction grows smoothly 
over a decade from zero to one, while in the latter it changes dramatically and almost instantaneously. Physically, 
what happens is that, at a certain temperature below the critical value, an enormous amount of hadronic bubbles 
are nucleated almost everywhere in the plasma, which grow explosively to transform all the plasma into hadrons. 
The picture is similar to what happen in standard cosmology [6] , in which the small supercooling phase transition is 
also very rapid with respect to the simple first order phase transition, at a temperature (T ~ 0.98rc)i which is very 
similar to the one we have obtained. Another remarkable feature that is present in the supercooling scenario is the 
important dependence of the transition temperature on the brane tension, about 10% difference between the lowest 
and the highest, which eventually leads to a different time at which the transition occurs. 
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V. DISCUSSIONS AND FINAL REMARKS 

In the context of brane-world scenarios, in the high density cosmological phase the Friedmann equation contains 
deviations to the AD case, which imposes fundamental phenonienological consequences on the cosmological evolution, 
and in particular on the cosmological phase transitions. In this work, we have analyze the evolution of the relevant 
physical quantities, namely, the energy density, temperature and scale factor before, during and after the phase 
transition. It is important to emphasize that in the early universe the energy density is extremely high, so that 
one can neglect the terms linearly proportional to the energy density with respect to the quadratic terms. This 
approximation has been considered in detail throughout this work. In the high density regime the Hubble function 
is proportional to the energy density of the cosmological matter, which drastically affects the cosmological dynamics 
of the universe. Moreover, in the early universe the dark radiation term, the projection of the Weyl tensor from the 
bulk, which appears in form of a radiation-like term in the field equations, may also play an important role. The 
magnitude of the brane world effects can be characterized by the numerical value of the brane tension A. In the limit 
A ^ oo we recover standard general relativity [5]. On the other hand, the study of the quark-hadron phase transition 
is also very important from an observational point of view, since the inhomogeneities generated at the QCD phase 
transition might have a noticeable effect on nucleosynthesis [9]. 

By fully taking into account the brane world effects we have found that the temperature evolution of the universe 
in the brane world scenario is different from the idealized standard FRW model. The temperature of the early 
universe in the quark phase is smaller in the brane world scenario, as one can see from Fig. 1, where the long dashed 
curve corresponds (approximately) to the general relativistic limit. Hence a small value of the brane tension would 
significantly reduce the temperature of the quark-gluon plasma, and accelerate the phase transition to the hadronic 
era. An increase of the dark radiation term for fixed bag constant B and brane tension A gives the same effect, 
as one can see from Fig. 2. Once the quark-hadron phase transition starts, the hadron fraction h is again strongly 
dependent on the brane tension. For small values of A, h(t) is much higher than in the standard general relativity, as 
one can see from Fig. 3. The increase of the dark radiation on the brane also strongly accelerates the formation of 
the hadronic phase (see Fig. 4), and decreases the time interval necessary for the transition. A small brane tension 
and a high energy density of the dark radiation also tend to reduce the temperature of the hadronic fluid. Of course, 
the temperature evolution also depends upon the relativistic degrees of freedom in the equation of state and upon 
the equation of state. In addition to this, by assuming that the phase transition may be described by an effective 
nucleation theory, we have also considered the case where the Universe evolved through a mixed phase with a small 
initial supercooling, and monotonically growing hadronic bubbles. It was shown that at a certain temperature below 
the critical value, an enormous amount of hadronic bubbles are nucleated, which grow explosively to transform all the 
plasma into hadrons, indicating that the small supercooling phase transition is very rapid with respect to the simple 
first order phase transition. 

Many details of the QCD phase transition are not yet conclusively understood. Even the order of transition is still 
a matter of debate. An advance in the understanding of the numerical values of the QCD coupling constants would 
be very helpful in obtaining accurate cosmological conclusions. Such an advance may also provide a powerful method 
for testing on a cosmological scale the theoretical predictions of the brane world models and the possible existence of 
the extra-dimensions. 
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